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It  is  conjectured  that  the  addition  of  a  small  compressive 
force  might  increase  bending  strength.  Extensive  tests  of 
western  hemlock  2  by  6’s  under  eccentric  axial  load  bear  this 
out.  Test  members  were  short  but  length  effect  is  studied 
analytically  by  simulating  long  members  with  a  finite  element 
computer  program.  A  small  sample  of  long  members  was 
tested  to  verify  the  computer  model.  Verification  was  not 
thorough  but  the  model  appears  to  do  very  well.  It  accurately 
reproduces  the  mean  but  underestimates  the  variance.  Further 
work  on  the  model  is  planned.  This  work  will  have  important 
applications  to  column  design  and  wood  truss  design. 
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Introduction 

Wood  members  under  combined  bending  and  compression  oc¬ 
cur  as  structural  elements  in  several  important  applications:  as 
the  top  chords  of  trusses,  as  wall  studs,  as  frame  members  in 
towers,  and  in  other  rigid  frame  structures.  While  strength  in 
bending  alone  and  compression  alone  have  been  extensively 
studied,  the  interaction  of  these  two  modes  of  failure  is 
poorly  understood,  particularly  for  dimension  lumber  of  com¬ 
mon  construction  grades.  There  is  reason  to  expect  that  a 
moderate  compressive  force  acting  in  concert  with  a  bending 
moment  could  raise  the  moment  capacity  by  shifting  the  mode 
of  failure  from  a  defect-controlled  brittle  rupture  to  a  more 
ductile  compression  failure  (/).’  Plastic  buckling  could 
become  a  failure  mode.  In  the  absence  of  good  data  on 
failure  under  combined  bending  and  compression,  present 
design  recommendations  are  quite  conservative;  this  conser¬ 
vatism  is  a  major  economic  factor  in  the  wood  roof  truss  in¬ 
dustry. 


Experimental 

i  Design  of  Experiments 

3>The  aim  of  this  study  was  to  define  a  failure  locus  in  the 
plane  of  bending  moment  (M)  versus  compressive  force  (C). 
Attention  was  restricted  to  that  portion  of  the  plane  where 
bending  stress  exceeds  compressive  stress  since  that  is  the 
region  of  greatest  importance  in  the  design  of  wood  frame 
structures  and  trussed  roofs. 


'  Maintained  at  Madison,  Wis.,  in  cooperation  with  the  University  of 
Wisconsin. 


In  order  to  limit  the  scope  of  the  study,  only  one  species 
(Western  Hemlock)  and  one  size  (nominal  2x6)  were  studied^ 
The  effects  of  length  and  grade  of  material  were  recognized  u  \ 
critical  but  rather  than  test  several  combinations  of  length  and 
grade  it  was  decided  to  simulate  these  effects  with  a  finite- 
element  computer  program.  In  that  way  only  short  pieces 
(single  “elements”)  need  be  tested  to  establish  a  data  pool 
from  which  a  computer  program  can  select  elements  and 
mathematically  simulate  the  behavior  of  longer  members.  The 
effect  of  material  grade  was  incorporated  by  measuring  a 
quality  index  for  each  element.  The  edgewise  bending  stiffness 
was  chosen  for  this  index  because  it  is  known  to  have  the  best 
correlation  with  bending  strength  of  any  nondestructive 
measurement.  In  addition  to  establishing  the  data  pool  of  ele¬ 
ment  properties,  a  few  longer  members  were  tested  to  verify 
the  finite  element  model’s  predictive  ability. 

Thus,  the  aims  of  the  experimental  program  were  twofold: 

1 .  To  test  short  members  under  combined  bending  and 
compression  at  various  M/C  ratios.  For  this  purpose  four 
large  groups  of  specimens  were  tested,  three  under  eccentric 
axial  load  at  different  eccentricities  and  one  in  pure  bending. 

2.  To  verify  the  computer-simulated  behavior  of  longer 
members.  For  this  (secondary)  purpose,  a  small  group  of 
8-foot  long  members  were  quality-rated  and  then  tested  under 
eccentric  axial  load. 

The  resulting  experimental  design  is  shown  in  table  1. 

Material 

Western  hemlock  of  10-foot  and  12-foot  lengths  was  pur¬ 
chased  from  the  Ebum  mill  of  the  Canadian  Forest  Products 


■  Italicized  numbers  in  parentheses  refer  to  literature  cited  at  end  of  this  report. 


1  These  had  been  pinned  to  lie  125  each,  but  terti  were  discarded 
whenever  failure  occurred  outside  the  gate  length. 


Company.  An  effort  was  made  toselect  members  with  a  Mo. 
2  grade  defect  in  the  gage  length  portions  of  two  5-foot 
members  to  be  cut  from  each  longer  length.  Since  western 
hemlock  and  white  fir  are  marketed  together,  each  piece  had 
to  be  identified  by  microscopic  examination  of  cell  structure 
and  the  white  fir  segregated  for  use  in  other  studies. 


Five-hundred  specimens  5  feet  long  were  cut  with  an  effort  to 
position  defects  in  the  center  18  inches  of  the  piece  whenever 
possible.  Thirty  members  11  feet  long  were  cut  from  pieces 
deemed  not  suitable  for  the  short  specimens.  Short  specimens 
were  quality  rated  by  testing  in  edgewise  bending  to  a  load  of 
300  pounds.  Two-point  load  was  applied  at  the  third  points  as 
shown  in  figure  1.  An  online  computer  took  digitized  readings 
of  load  and  head  movement  at  two  points  and  computed  an 
“E-rating”  value,  E,,  from  the  formula 
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Figure  2.— Loading  used  for  E-rating  specimens  in  group 
five.  The  member  was  rated  five  times,  one 
for  each  19-inch  element.  The  position  shown 
here  is  for  rating  element  2. 


E,=  5(Pj-P,)L’ 

27{H,-H,)bh* 

where  P^load,  H«=head  movement,  L=»=54  indies,  b= width, 
and  h= depth.  The  11-foot  specimens  were  each  quality  rated 
at  five  different  positions  by  a  similar  test  in  edgewise  bending 
under  two-point  loading  as  shown  in  figure  2.  Between 
readings  the  member  was  indexed  through  the  test  machine  in 
19-inch  increments.  Again  a  computer  used  head  movement 
and  machine  load  to  compute  E,.  for  each  of  five  “elements" 
in  the  central  8-foot  gage  length  of  these  11-foot  members. 


1  Maan  value  (coefficient  of  variation). 

*  Dorn  am  apply  (Group  4  was  tested  la  bmdiag). 
’  Not  obtained, 

HWwH  wt  cconuiuiy*  >*v*a  h  ■  mmwiq 

•  E.  was  of  the  earner  dement  In  Group  S. 


from  the  Has  of  action  of  the  axial  load. 


Figure  3.— Loading  fixture  for  applying  eccentric  axial 
load:  A)  diagram,  B)  photo  of  side  view,  and 
C)  photo  looking  into  empty  boot.  Rope 
restrains  boot  from  falling  out  when  specimen 


The  short  members  were  sorted  in  order  of  increasing  Er.  The 
first  four  were  randomly  assigned  to  groups  1  through  4.  The 
next  four  were  then  randomly  assigned,  and  the  next  four, 
etc.,  until  four  matched  groups  were  produced  with  nearly 
identical  distributions  of  Er.  Table  2  shows  how  closely  the  Er 
distributions  of  groups  one  to  four  agreed  Means  ranged 
from  1.03  to  1.03  million  pounds  per  square  inch  and  coeffi¬ 
cients  of  variation  ranged  from  0.137  to  0.211. 


A!=  jfy,  +y2)  I 

1  =  2(y,-y,)  I 

R  a(2!  +y,  +y2) 

where  A i  =  axial  elongation  (contraction  if  negative),  R 
-radius  of  curvature,  t  =  gage  length  (18  in.),  and  a  =arm 
length  (14  in.). 


Test  Methods 
Groups  One  to  Three 

Eccentric  axial  load  was  applied  through  grips  in  the  form  of 
steel  boots  18  inches  long.  The  top  of  the  boot  bore  against 
the  end  of  the  specimen.  The  narrow  sides  of  the  boot  were 
made  snug  against  the  specimen  by  means  of  adjusting  screws 
shown  in  figure  3.  The  bottom  of  the  boot  had  flanges  which 
received  the  test  machine  load  by  way  of  a  yoke  which 
straddled  the  boot.  This  put  the  pivot  point  as  close  to  the 
gage  length  as  possible  in  order  to  minimize  the  growth  of 
eccentricity  due  to  specimen  deformation  during  a  test.  Since 
inservicc  wall  studs  and  roof  truss  chords  are  restrained 
against  lateral  buckling  by  sheathing  materials,  the  test 
specimens  were  laterally  supported  at  midspan  by  fixed  braces 
on  each  side  of  the  specimen.  Initially  there  was  an  1 /8-inch 
clearance  between  the  specimen  and  the  lateral  supports  and 
the  support  braces  were  surfaced  with  a  hard  polished  plastic 
having  a  low  coefficient  of  friction.  In  this  manner  the 
specimens  were  constrained  to  deflect  only  in  the  direction  of 
greatest  stiffness. 

Deformation  measurements  were  taken  by  two  linear  variable 
differential  transformers  (LVDT’s)  at  the  ends  of  arms  at¬ 
tached  to  the  wide  faces  of  the  specimen,  as  shown  in  figure 
4.  The  arms  were  clamped  to  the  specimen  with  an  18-inch 
gage  length  between  them.  The  two  deflection  readings  y,  and 
yi  are  related  to  axial  deflection  and  curvature  as  follows: 
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Figure  4.  —Deflectometer  arrangement.  Rigid  arms  were 
clamped  to  specimen  and  linear  variable  dif¬ 
ferential  transformers  (LVDT’s)  were  mounted 
between  them  at  7  inches  on  either  side  of 
center  line. 


M  140642 


Figure  5. — Exaggerated  specimen  deformation  showing 
curvature  and  center  deflection,  6 


The  moment  applied  to  the  specimen  by  the  eccentric  axial 

te-U  M-P4  (4) 

where  6  equals  lateral  deflection  of  specimen  centroidal  axis 
measured  from  line  of  action  of  P.  During  test,  the  boots 
rotate  through  an  angle  0 


where  S  equals  arc  length  between  boots,  which  is  24  inches 
(see  figure  5).  Thus  A  can  be  computed  from 

d=e  cos-®  +  b  sin|+R(l-cosj)  (6) 

where  e- eccentricity  and  b=4  inches = distance  from  pivot 
point  to  edge  of  boot. 

Data  Recording 

Load  P  and  deflections  y,  and  y2  were  recorded  automatically 
on  punched  paper  tape  with  the  aid  of  an  online  computer 
and  teletype  machine  (see  figure  6).  An  X-Y  plotter  recorded 
y,  and  y,  versus  P.  The  paper  tapes  were  later  read  by 
another  computer  and  reduced  to  R,  P,  and  6  and  stored  on' 
magnetic  tape.  These  magnetic  tapes  and  the  X-Y  plots  form 
a  permanent  record  of  the  raw  test  data. 


*1 
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Figure  6.— Automatic  data  recording  was  accomplished 
with  (a)  an  X-  Y  plotter,  fb)  an  online  compu¬ 
ter,  and  (c)  a  teletype  machine  with  papa’  tape 
punch. 
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Data  Prow li 

Each  test  in  groups  one  to  four  was  processed  to  yield  a  plot 
of  moment  versus  curvature.  From  these  plots  four  properties 
were  extracted:  stiffness,  strength,  ductility,  and  knee-shape. 
To  do  this,  the  Ramberg-Osgood  function  was  fitted  to  the 
plots  of  M  versus 


where  c = half  depth = 2.73  inches,  El  =  flexural  rigidity, 
pounds-inchesj  and  K  and  n  are  fitted  parameters. 

To  minimize  rounding  error  in  the  curve  fitting  computations 
it  is  better  to  write  equation  (7)  in  the  form 

^ = x + exp  [ln(y„) + n  ln(|-) )  (8) 

where  x  *  gjc,  xu  ■  x  at  ultimate  load,  and  y„  *  at 

ultimate  load.  This  reduces  to  the  Ramberg-Osgood  form  (eq. 
(7) )  if  K  *  K  is  not  a  convenient  parameter  since  it 

can  easily  become  so  large  that  it  overflows  the  floating  point 
accumulator  of  a  digital  computer.  Fitting  the  function  given 
by  equation  (8),  the  four  desired  properties  can  easily  be 
identified  with  the  following  fitted  parameters: 

El  =  stiffness,  xu  =  strength,  yu  =  ductility,  and  n  =  knee-shape. 

Figure  7  shows  a  typical  fitted  curve  and  the  geometric  mean¬ 
ing  of  the  four  properties.  Hereafter  this  moment-curvature 
characteristic  will  be  referred  to  as  the  “property”  curve  since 
it  is  uniquely  specified  by  the  four  relevant  properties.  Note 
that  the  tail  of  the  test  record  beyond  maximum  load  is  ig¬ 
nored  in  fitting  the  Ramberg-Osgood  function,  because  this  is 
a  test  machine  artifact  and  not  relevant  to  the  dead-load  con¬ 
ditions  for  which  members  are  designed. 
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Figure  7.— Typical  plot  of  moment  versus  curvature, 

from  which  four  element  properties  can  be  de¬ 
rived:  stiffness  El,  strength  Mv  ductility  yw 
and  shape  n.  Larger  n  implies  a  sharper  knee. 


Minor  Testa 

After  testing  the  strength  specimens  to  destruction,  small 
coupons  approximately  1.3  by  1.3  by  3.3  inches  were  cut  from 
each  and  used  to  measure  specific  gravity  and  moisture  con¬ 
tent.  The  results  are  summarized  for  each  group  in  table  2 
which  shows  means  followed  by  coefficient  of  variation  in 


parentheses.  It  is  seen  that  these  properties  did  not  vary  ap¬ 
preciably  between  groups.  Therefore  no  adjustments  woe 
made  to  the  data  to  correct  for  these  variations. 

Results  of  Experiments 

Short  Specimens — Groups  One  to  Four 

The  results  of  the  short  member  tests  were  four  properties  of 
each  specimen  obtained  by  fitting  a  Ramberg-Osgood  func¬ 
tion  as  described  under  Data  Processing.  There  were  four 
groups  of  specimens  corresponding  to  four  eccentricities  as 
shown  in  table  1.  Within  each  group  a  normal  probability 
curve  was  fitted  to  each  of  the  four  properties,  stiffness  El, 
strength  xu,  ductility  yu,  and  knee-shape  n,  and  to  the  quality 
index,  Er.  In  some  cases  it  was  found  that  the  logarithm  of 
the  property  was  more  nearly  normal  than  the  property  itself. 
Specifically,  the  following  properties  were  fitted  with  a  multi¬ 
variate  normal  probability  curve: 

Stiffness— El,  Strength— In  Mu,  Ductility— In  yu.  Knee- 
shape— In  n,  Quality— Er,  where  Mu  is  the  moment  at  max¬ 
imum  load  and  E,.  is  the  E-rating  index  of  quality  obtained 
from  a  preliminary  nondestructive  test  as  described  under 
Specimen  Preparation.  The  results  of  these  curve  fittings  are 
summarized  in  table  3  which  gives  the  means  and  the  matrix 
of  covariance  for  groups  one  to  four. 

Since  the  covariance  matrix  is  diagonally  symmetric,  it  con¬ 
tains  only  IS  independent  elements.  These  plus  the  3  means 
form  20  independent  statistics  for  each  group.  These  shall  be 
denoted  as  S,  through  S»  defined  such  that 


El 

s, 

ln(Mu) 

s. 

ln(yu) 

= 

s, 

ln(n) 

s. 

3  . 

s, 

S.  S, 


Covariance  matrix  = 
(diagonally  symmetric) 


•  Si, 

•  • 

•  • 


•  • 


S.  S*  Si« 
S|2  S„  S|4 
Sis  Si*  S|7 

•  Sib  Sib 

•  •  Sjo 


For  example,  St  is  the  covariance  of  El  and  In  yu.  Together, 
the  20  statistics  for  any  one  group  uniquely  specify  a  multi¬ 
variate  (specifically  five-variate)  normal  distribution  that 
characterizes  that  group. 


It  should  be  noted  that  the  characteristic  multivariate  distribu¬ 
tion  depends  on  e,  the  eccentricity  of  the  group.  All  20 
statistics  (with  the  exception  of  Sj  and  S»,  the  mean  and 
variance  of  Er)  are  functions  of  eccentricity  e.  Figure  8 
shows,  as  an  example,  how  the  mean  of  the  strength 
parameter  In  Mu  varies  with  eccentricity.  For  convenience  in 
treating  group  four  (pure  bending,  e  =  °°)  the  statistics  are 
plotted  against  the  reciprocal  of  e  and  group  four  is  at  the 
origin.  Note  that  the  expected  results  were  obtained:  pure 
bending  does  not  have  the  highest  moment  capacity.  It  is  in¬ 
teresting  to  plot  this  in  a  plane  of  moment  versus  axial  force. 
Figure  9  shows  such  a  plot.  It  is  seen  that  the  addition  of  a 
small  amount  of  axial  compression  does  indeed  increase  the 
moment  capacity  for  members  tested  over  a  short  gage  length. 
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Figure  8.— Variation  of  input  distribution  with  ecentri- 
city.  The  mean  of  In  Mu  is  shown,  along  with 
a  jilted  parabola. 


Long  Specimens— Group  Five 

It  is  tempting  at  this  point  to  simply  use  the  data  of  groups 
one  to  four  as  a  basis  for  design  of  members  under  combined 
bending  and  compression.  However,  the  effect  of  length  (in¬ 
cluding  the  possibility  of  plastic  buckling)  is  critical  and  is  not 
contained  in  figure  9.  Group  S  shows  this.  It  had  the  same  ec¬ 
centricity  as  group  one  but  different  length  and,  as  table  2 
shows,  its  mean  ultimate  load  was  only  about  half  of  that  of 
group  one.  The  Er  distributions  of  groups  five  and  one  are 
not  quite  as  closely  matched  as  those  of  groups  one  to  four 
but  are  still  very  close.  This  comparison  clearly  shows  that 
there  is  a  profound  length  effect.  It  is  due  to  plastic  instability 
in  the  presence  of  material  asymmetries  around  defects  and 
variable  stiffness  along  the  length.  Such  an  effect  is  difficult 
to  analyze  mathematically  in  closed  form.  A  nonlinear, 
numerical  analysis  is  done  in  the  next  section,  where  the  Er 
variation  along  the  length  of  the  members  in  group  five  is 
used  as  an  input  to  a  non-standard  finite  element  simulation 
of  the  group  five  specimens. 

Monte  Carlo  Simulation  of  Length  Effect 

Group  five  specimens  were  simulated  with  the  special  finite 
element  analysis  given  in  the  appendix.  In  that  analysis  the 
nonlinear  Ramberg-Osgood  function  given  in  equation  (7)  is 
used  to  represent  the  relationship  between  moment  and  curva¬ 
ture.  This  equation  is  exactly  solved  by  iteration  for  each 
given  set  of  four  Ramberg-Osgood  properties:  El,  xu,  yu,  and 
n.  These  properties  are  assumed  to  be  piecewise  constant 
functions  of  x  with  five  “pieces”  or  “finite  elements”  cover¬ 
ing  the  range  of  x.  Deflections  and  slopes  are  matched  at  the 
nodes  between  elements.  This  is  far  more  exact  than  a  stan¬ 
dard  finite  element  approach  which  would  approximate 
deflection  shapes  over  each  element  and  would  account  for 
geometric  nonlinearity  in  an  approximate  manner  via  the  so- 
called  stability  matrix.  Since  the  number  of  elements  could 
not  be  made  large  here,  it  was  decided  to  use  a  more  exact 
non-standard  finite  dement  method.  Five  19-inch  elements 


Figure  9.— Results  of  short  member  tests:  failure  moment 
versus  failure  load  (both  measured  at  point  of 
maximum  toad!  for  groups  one  to  four 
(KIPS  =■=  thousand  lbs). 

add  up  to  a  95-inch-long  member  —  essentially  8  feet.  The 
results  of  the  short  member  tests,  which  had  a  gage  length  of 
18  inches,  were  used  in  a  five-element  model  to  simulate  the 
behavior  of  the  group  five  specimens.  The  error  of  using  data 
obtained  over  18  inches  in  a  19-inch  element  is  negligible. 

Input  Parameters 

In  the  Monte  Carlo  technique,  input  parameters  are  randomly 
selected  from  a  known  distribution  and  a  test  is  simulated. 
This  is  repeated  many  times  until  enough  output  strength 
values  have  been  generated  to  give  an  idea  of  the  output 
distribution. 

In  this  case  the  input  was  of  two  kinds:  (1)  the  distribution  of 
element  properties  and  (2)  a  quality  description  of  the  long 
member  to  be  simulated.  Input  (1)  consisted  of  a  five-variate 
normal  distribution  in  the  variable  El,  ln(M  ),  ln(yu),  ln(n), 
and  Er  This  was  specified  by  a  five-vector  of  means  and 
a  five  by  five  symmetric  covariance  matrix,  that  is,  by  20 
independent  statistics.  See  equations  (9)  and  (10).  With  the 
exception  of  Ef  and  its  variance,  each  of  these  statistics  is  a 
function  of  the  eccentricity  of  the  element  as  can  be  seen  by 
examining  table  3  and  figure  8.  Hence  each  was  specified  by 
the  three  coefficients  of  a  parabola 

Si  =  ai(i)i  +  bi(i)  +  ci  (11) 

Here  aj,  bj,  and  Cj  are  the  fitted  coefficients,  Sj  is  a  statistic 
fitted  by  least  squares  to  the  data  in  table  3,  and  e  is  the  ec¬ 
centricity.  These  coefficients  are  shown  in  table  4.  The  other 
input  (2)  consisted  of  the  E,.  values  of  the  five  elements  of  the 
long  member  to  be  simulated.  These  had  been  obtained 
nondestructive^  for  the  specimens  in  group  five  before  testing 
them  in  eccentric  axial  compression. 
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Random  Selection  of  Elements 

The  Monte  Carlo  technique  proceeded  as  follows: 

For  each  element  with  a  given  measured  value  of  Er, 

1.  Generate  four  independent  random  normal  variates  with 
zero  mean  and  unit  variance: 


This  vector  will  be  used  in  step  4.  It  identifies  the  randomly- 
chosen  element. 

2.  Using  the  current  eccentricity  e  of  the  element,  calculate 
the  parameters  S,  through  S»  of  the  element  property 
distribution  at  this  eccentricity: 


ln(Mu) 

ln(yu)  =  Normal  with  mean  and  covariance  X 
ln(n) 

Er  J 

where 


Table  3. — Seminary  of  property  statistics 


Statistic 

Group 

1 

2 

3 

4 

MEANS  (see  eq.  (9» 

Si  =  mean(KI) 

0.444761+ OS 

0.400781+08 

0.391585  +  08 

0.448324+08 

S,  =  meant  In  Mu) 

.104730+02 

.105860  +02 

.107590+02 

.106890  +  02 

S,  =  meant  In  yu) 

-.687770+01 

-.672110+01 

-.688740  +  01 

-.787720  +  01 

S.  =  meantln  n) 

.205100  +  01 

.219090+01 

.209800  +  01 

.165650+01 

Ss  =  mean(Kr) 

.104183  +  07 

.104678+07 

.102889  +  07 

.104045  +  07 

COVARIANCES  (see  eq.  (10)) 

S.  *var(EI) 

.133746+15 

.106474+15 

.938794  +  14 

.107557  +  15 

S,  =  cov(El,  In  Mv) 

.138715  +  07 

.153167+07 

.227157  +  07 

.326864  +  07 

St  -co»(IE,  la  yn) 

-.405308  +  07 

-.179065  +  07 

-.176640  +  07 

.378817  +  06 

S,  =  cov(F.I,  la  a) 

.115347+07 

.217908  +  07 

.256852  +  07 

.191873+07 

Sit  -covtn.EjJ 

.106694  +  13 

.126682  +  13 

.153258  +  13 

.141881 + 13 

Si,  -vartlaMJ 

0.632824-01 

0.578820  -  01 

0.116250  +  00 

0.192532+00 

S„  =cov(ln  Mu,  la  yj 

.473633  -  01 

.309349  -  02 

.291294  -  01 

.772557  -  01 

Si,  =  cov(ln  Mu>  In  a) 

.243918  -  01 

.456674  -  01 

.116520  +  00 

.885921-01 

S,t  =  cov(tn  Mu,  E^ 

.254907  +  05 

.274862  +  05 

.470426  +  05 

.444307  +  05 

S„  =  vartln  y^) 

.685311+00 

.451765  +  00 

.106489  +  01 

.995538+00 

S„  =cov(ln  yu,  la  a) 

.420844  -  01 

-.145156  -  01 

.153682  +  00 

.703405  -  01 

S„  =cov(lay  n,  E^ 

-.298365  +  05 

-.225922  +  05 

-.351988  +  05 

.915227  +  04 

S,t  =  var(la  a) 

.217163+00 

.328039  +  00 

.547928  +  00 

.386828  +  00 

S.t  »cav(ln  a,  Ep 

.957648  +  04 

.174911+05 

.558619  +  05 

.337252  +  05 

Sw  -vartEjJ 

.352101  +  11 

.314806+11 

.469594  +  11 

.268380+11 

and 


Table  4. — Coefficients  of  fitted 


(seeeq.  (11» 


s. 

s, 

s. 

S., 

s, 

Su 

Su 

S„ 

I-  s, 

s„ 

Su 

Su 

s. 

s„ 

S,. 

Si» 

^SlO 

Su 

s„ 

Si* 

where  Sj,i  =  1,20  are  given  by  equation  (11)  as  functions  of  e. 
(Note  that  S,,  and  S„  are  independent  of  e,  so  the  given  Ej  is 
always  constant.) 

3.  With  this  p  and  I  and  the  given  Er,  calculate  a  condi¬ 
tional  four-variate  property  distribution: 


=  Normal  with  mean  £iand  covariance  £ 


+  hzh 

S]0 


Si. 

S,*S,4 

S|tS|7 

s,.s 

S|oS|4 

S|4 

SI4S17 

SuS 

*  S|,S„ 

S|?S|4 

S?7 

S„S 

Si*St* 

S|*S|4 

S|*S,7 

Si. 

4.  Calculate  the  four  element  properties 


El 

ln(Mu) 

ln(yu) 

ln(n) 


Z  +  M 


From  these  four  values  the  property  curve,  equation  (8),  of 
the  element  can  be  determined. 


Element  Identity 

As  load  increases  the  eccentricity  of  any  given  element  in¬ 
creases  in  response  to  the  increasing  deflections  of  the 
member.  And  since  it  was  found  that  characteristic  element 
behavior  is  different  at  different  eccentricities  it  was  decided 
to  let  the  characteristic  property  distribution  of  the  element 
evolve  as  a  function  of  where  e  is  the  current  eccentricity 
of  that  element.  Thus,  element  properties  are  not  constant 
and  cannot  be  said  to  identify  the  element.  However,  the  z 
vector  can  be  used  to  identify  the  element,  since  it  specifies 
the  number  of  standard  deviations  away  from  the  mean  that 
the  element  lies  within  its  distribution.  If  this  is  held  constant, 
the  element  keeps  its  rank  order  relative  to  other  elements  in 
the  distribution.  This  is  like  assuming  that  the  rank  order  of 
the  short  specimens  in  any  group  would  have  been  the  same  if 
they  had  been  tested  at  any  other  eccentricity.  (This  assump¬ 
tion  is  a  useful  expedient,  but  it  may  not  be  strictly  te,  since 
the  relative  severity  of  defects  may  be  influenced  to  x- 
tent  by  the  eccentricity  of  the  load.) 


“  I  "1 

f.'A 


i 

“i 

bi 

ci 

1 

0.233530  +08 

-0.232150+00 

0.445570  +  08 

2 

-.384481  -f  00 

.135020+00 

.107071+02 

3 

-.274490  +  01 

.365600  +  01 

-.702813+01 

4 

-.139041+01 

.175730+01 

.167184+01 

s 

0.0 

0.0 

.103719+07 

4 

.922241  +  14 

-.646011  +  14 

.106760+15 

7 

.220446+07 

-.400230  +  07 

326550+07 

• 

-.131435+06 

-.374509+07 

.702772  +  05 

9 

-380925+07 

.290017+07 

.194974+07 

1* 

-.666018+12 

.266009  +  12 

.144551  +  13 

11 

.207180+00 

-.330225  +00 

.193529+00 

12 

.222322  +00 

-.255763  +  00 

.792267  -  01 

13 

-.900562-01 

.177533-01 

.975369-01 

14 

-.605011+04 

-.170433+05 

.467426+05 

15 

322511+00 

-.994909  +  00 

.108568+01 

M 

-.261666-01 

-340651-01 

.996603-01 

17 

.052610  +  05 

-.114556+06 

376651+04 

11 

-303959  +00 

358054+00 

.417622+00 

19 

-.613004+05 

.271180+05 

.393404  +  05 

20 

0.0 

0.0 

304231  +  11 

M14SS73 

Figure  10.— Block  diagram  of  computing  method. 
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Computing  Method 

Figure  10  shows  a  block  diagram  of  the  computing  method. 
The  inner  loop  is  necessitated  by  the  nonlinearity  of  the 
moment-curvature  relationships  of  the  elements.  The  moment- 
curvature  relationship,  or  “property  curve,”  is  represented 
here  as  a  Ramberg-Osgood  function  as  discussed  under  Data 
Processing.  As  the  deflections  increase,  the  property  curves 
must  shift  in  response  to  the  increasing  eccentricities  of  the 
elements.  This  necessitates  a  middle  loop  of  iterations  until 
the  deflections  of  the  nonlinear  finite  element  solution  are 
congruent  with  the  eccentricities  used  to  derive  the  properties 
of  each  element.  Finally,  there  is  an  outer  loop  in  which  the 
load  is  gradually  incremented  until  at  least  one  element 
reaches  failure.  If  these  increments  are  not  too  large,  conver¬ 
gence  of  the  inner  loops  will  be  very  rapid.  A  complete  flow 
diagram  of  the  subprogram  SIMTST  used  to  simulate  an  ec¬ 
centrically  loaded  member  is  shown  in  figure  11.  Computing 
time  to  simulate  one  destructive  load  test  of  a  five-element 
member  was  about  4  seconds  on  a  Univac  1 1 10  computer. 

Results  of  Monte  Carlo  Simulations 


Figure  12. — Simulated  failure  load  of  8-foot-long  mem¬ 
bers  versus  experimental  failure  toads  of 
group  five. 


Verification  of  the  Model 

Figure  12  shows  a  comparison  of  model  and  experimental 
results.  Group  five  consisted  of  a  small  sample  of  just  29 
members  of  8-foot  gage  length  under  1-inch  eccentric  axial 
load.  This  group  was  included  in  the  study  in  order  to 


Figure  11.— Flow  diagram  for  program  SIMTST,  the 
finite  element  simulation  of  a  long  member 


test. 


roughly  verify  the  computer  finite  element  model.  Each  test 
member  had  been  E-rated  at  five  locations  within  the  gage 
length,  and  the  computer  used  these  ratings  as  input  specifica¬ 
tions  for  a  five-element  simulation.  Each  member  was 
simulated  30  times  with  element  properties  chosen  randomly 
from  conditional  distributions,  given  the  input  E-ratings.  The 
simulations  are  shown  in  figure  12  as  vertical  lines  represent¬ 
ing  the  mean  plus  or  minus  one  standard  deviation.  The  45  0 
lines  show  where  the  data  should  fall  for  perfect  verification. 
It  is  seen  that  the  mean  of  group  five  is  accurately  reproduced 
but  that  the  variance  is  underestimated.  That  is,  the  trend  of 
the  data  is  not  a  45  °  line  as  it  should  be.  This  suggests  that 
there  may  be  a  better  index  of  quality  than  the  edgewise 
E-rating.  Perhaps  flatwise  E  would  prove  to  be  more  sen¬ 
sitive.  The  correlation  coefficient  between  Er  and  In  Mu  was 
on|y  [b.54]  IT 

correlation  coefficient  =  for  group  ^  (18) 

0.60  4 


in  this  study,  which  is  not  very  great.  The  correlation  between 
the  E,.  of  the  center  element  and  the  ultimate  load  for  group 
five  was  also  small,  namely  0.68.  But  the  fact  remains  that  the 
correlation  between  strength  and  nondestructive  quality  in¬ 
dices  is  never  very  great.  This  is  the  essential  problem  of 
lumber  grading  translated  here  into  element  grading.  Never¬ 
theless  it  may  be  possible  to  force  the  computer  model  to  do 
better  in  the  aggregate  even  though  it  remains  a  somewhat 
poor  predictor  for  individual  members.  This  is  worthy  of  fur¬ 
ther  study. 


The  Effect  of  Length 

The  effect  of  length  has  been  discussed  briefly  above  in  the 
comparison  of  experimental  results  for  groups  one  and  five. 

A  more  thorough  investigation  of  the  effect  of  length  can  be 
made  with  the  use  of  the  computer  model.  Figure  13  shows 
how  the  29  specimens  in  group  five  would  have  behaved  if 
they  had  been  tested  at  other  eccentricities  or  in  pure  bending. 
Each  specimen  is  represented  by  the  average  of  30  simulations 
and  mean  and  standard  deviation  of  the  sample  of  29  is 
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Figure  13. — Computer  simulations  of  8-foot-long  mem¬ 
bers  under  eccentric  axial  load. 


£>,  IO*psi 
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Figure  14. — Comparisons  of  Er  distributions  of  group 
five  and  groups  one  to  four. 


plotted  in  the  figure.  Deflections  due  to  load  are  ignored  in 
plotting  moment  versus  axial  load,  that  is: 


M  =  Pe 

is  used  rather  than 

M  =  Pd 


(19) 

(20) 


even  though  the  deflection  6  might  be  quite  large.  This  was 
done  for  the  sake  of  facilitating  application  to  design.  The 
design  engineer  knows  the  eccentricity  e  but  cannot  accurately 
estimate  the  late,  at  deflection  A  especially  when  plastic 
buckling  is  a  possibility. 


It  is  seen  in  figure  13  that  the  mean  moment  capacity  does 
not  improve  with  the  addition  of  a  small  amount  of  axial 
compressive  force  as  was  found  for  groups  one  to  four. 
However,  the  low  tail  of  the  distribution  does  still  exhibit  this 
effect,  though  it  is  far  less  pronounced  than  it  was  for  the 
short  specimens  (compare  fig.  9).  Before  this  effect  can  be 
counted  upon  in  design,  other  loadings  and  different  species 
and  sizes  should  be  investigated.  Also,  the  behavior  of  the 
computer  model  in  the  low  tail  may  not  be  reliable  since  the 
model  is  known  to  underestimate  variability.  Of  course,  for 
very  long  members,  failure  by  buckling  will  appear  at  low 
compressive  loads.  Even  at  8  feet  the  influence  of  some  in¬ 
elastic  buckling  is  already  evident,  as  can  be  seen  by  compar¬ 
ing  figures  9  and  13.  In  this  comparison  it  is  important  that 
the  grade  of  the  small  population  in  group  five  should  resem 
ble  the  grade  of  the  populations  in  groups  one  to  four.  Figure 
14  shows  a  comparison  of  the  Er  distribution  of  group  five  to 
the  common  Er  distribution  of  groups  one  to  four.  It  is  seen 
that  the  resemblance  is  quite  close  although  no  special  attempt 
was  made  to  select  group  five  members  for  this  match.  Also, 
the  comparison  of  figures  9  and  13  is  a  comparison  of 
simulated  members  with  real  test  specimens,  so  modeling  error 
is  present.  The  model  underestimates  variability.  Thus  the 
comparison  should  be  restricted  to  the  means.  Nevertheless, 
the  comparison  strongly  indicates  that  the  effect  of  length  is  a 
profound  one.  Three  things  contribute  to  this:  inelastic 


buckling,  size  effect,  and  the  interaction  of  bending  and  com¬ 
pression  failure. 

Summary 

The  interaction  of  bending  and  compression  strengths  of 
western  hemlock  2  by  6  lumber  was  investigated.  Over  400 
short  specimens  (with  gage  length  =18  in.)  were  tested  in  four 
groups  with  different  ratios  of  moment  and  compression. 
Some  interaction  was  found  that  decreases  with  increasing 
length,  as  shown  by  the  results  of  a  finite  element.simulation 
also  reported  here.  This  finite  element  model  was  able  to 
reproduce  mean  results  of  a  small  sample  of  8-foot-long 
specimens  but  did  not  show  as  large  a  variation  as  the  experi¬ 
mental  data.  The  model  could  be  greatly  improved  if  a  better 
nondestructive  index  of  quality  could  be  found  which  cor¬ 
related  better  with  strength  than  the  edgewise  bending  stiff¬ 
ness  used  here.  Perhaps  a  combination  of  two  indices  should 
be  used. 

Data  should  be  obtained  for  other  sizes  and  species.  The 
results  could  be  used  by  researchers  studying  the  reliability  of 
wood  walls  and  roof-truss  systems.  The  data  presented  here 
are  sufficient  to  show  that  a  simple  linear  interaction  failure 
criterion  is  conservative — the  degree  of  conservatism  being 
greater  for  shorter  members. 
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Appendix:  Finite  Element  Analysis 


Combined  Bending  and  Compression 

Consider  a  member  of  length  L  under  eccentric  axial  load  and  divide  the  member  into  m  equal 
dements  of  length  h  as  shown  in  figure  Al. 


L  =  mh 

(Al) 

Isolate  the  nth  element.  From  elementary  beam  theory  the  moment-curvature  relationship  is 

<E1>ny;  =  -pyn 

(A2) 

where  yn  =  lateral  deflection  in  the  range  (n-1)  h  <  x  <  nh  measured  from  the  line  of  action  of  P 
and  primes  denote  differentiation  with  respect  to  x.  Let  the  nodal  displacements  be  q„: 

y  =  qn,  at  x  =  nh,  n  =  0,m 

(A3) 

The  solution  of  differential 

equation  (A2)  which  satisfies  boundary  conditions  (eq.  (A3))  is 

V  s 

qncos(n  -  l)Ajjh  -  qn  _  jCOS(nA„h) 

sinA nx 

7n 

sin^h 

+ 

["  -q„sin(n-  OAjjh  +  q„_  jSinfiiAjjh) 

cosAjjX 

(A4) 

i  sinAjjh 

where 

/  P 

*n“  V  (EI)n 

Match  slopes  at  nth  node: 

y’n  =yn  +  l  =  0atx  =  nh 

(A6) 

Substituting  equation  (A4)  into  (A6)  and  using  some  trigonometric  identities  this  becomes 


*n 

V°SV>  ( 

An  +  jC°sAn  +  1h 

*n  +  l 

s^h 

qn-l  + 

sinAjjh 

sinAn+,h 

q„  +  - 

sinAn+,h 

Equation  (A7)  yields  a  banded  matrix  of  m-1  equations  for  m  + 1  unknown  q0  to  qm.  But,  of 
course,  two  of  the  nodal  displacements  are  given: 


do  -  =  e 

where  t  is  the  given  eccentricity.  Thus  all  nodal  displacements  can  be  obtained  by  solving  an 
(m-l)th  order  linear  system  if  the  A„  are  given.  For  small  loads,  the  (EI)n  are  constants  indepen¬ 
dent  of  P.  For  larger  P  one  must  obtain  (EI)n  from  the  secant  to  the  point  on  the  property  curve 
corresponding  to  the  load.  To  find  this  point  requires  some  iteration,  since  the  deflections  de¬ 
pend  on  the  (EI)n,  the  (El),,  depend  on  the  moments,  and  the  moments  depend  on  the  deflec¬ 
tions.  A  few  iterations  around  this  loop  will  quickly  converge,  however,  especially  if  one  starts 
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Figure  A 1.— Fret-body  diagram  of  member  under  eccen¬ 
tric  axial  load.  Centroida I  deflection  y  is 
measured  from  line  of  action  of  load  P. 
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Figure  A3.  -Free-body  diagram  of  member  under  pure 
bending. 


Figure  A2.  —Iteration  loop  used  in  finite  element  simula¬ 
tion  program. 


with  a  good  initial  guess.  Hence  for  large  loads  one  must  begin  with  a  load  small  enough  to  en¬ 
sure  linear  behavior  and  increase  by  small  increments,  iterating  to  convergence  after  each  incre¬ 
ment,  and  using  the  solution  from  the  last  load  as  the  initial  guess  for  the  (EI)„  after  each  incre¬ 
ment.  This  procedure  continues  until  one  of  the  elements  reaches  its  failing  load.  Figure  A2 
shows  the  iteration  loop. 

This  method  is  capable  of  accurately  modeling  the  phenomenon  of  plastic  buckling  when  all 
elements  are  sufficiently  ductile.  In  that  case  the  deflections  grow  rapidly  just  before  failure. 

Pore  Bending 

In  the  case  of  pure  bending  the  finite  element  analysis  is  greatly  simplified  by  the  absence  of  any 
buckling  phenomenon.  Cbnsider  a  member  of  length  L  under  pure  bending  as  shown  in  figure 
A3.  Again  cut  the  member  into  m  elements  of  length  h  so  that 

L  =  mh  (A9) 

Isolate  the  nth  element.  The  moment-curvature  relationship  is 

mnyn  =  M  (A10) 

where  yn  =  lateral  deflection  in  the  range  (n  -  l)h  <  x  <  nh.  Let  nodal  displacements  be  qn: 

y  =qn  at  x  =  nh,  n  =  O.m  (All) 

The  solution  of  equation  (A  10)  satisfying  conditions  (All)  is 

yn  “  f">  “  Mnh(n-i)x  +  p„(^)n(n-l) 

+  q„(g)  -  (r-  -  l)qn  -  qn  _  ,(*)  +  nq„  _  ,  (A12) 

where 

■  M 

(EI)n  (A13) 


Match  slopes  at  nth  node: 


yn  -  yn  + 1  -  °-x  *  nh 


(A14) 
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Substituting  (A12)  into  (A14)  and  simplifying  yields 


%-l  -  2%  +  %-t-l  "  rW  +^+1)  <A,5> 

This  is  a  finite  difference  equation  whose  solution  can  be  obtained  in  closed  form.  The  homoge¬ 
neous  solution  is 

qn  =  a  +  pn  (A16) 

where  a  and  jJ  are  arbitrary  constants.  Now  let  a  and  fi  be  functions  of  n: 

%  =  tfn  +  V  <A17> 

and  use  the  method  of  variation  of  parameters  ( 2 )  to  find 

n  - 1  2 

«„  =  -I  0  +  D  jOtj+Pj  +  f)  +  A  (A18) 


Pa  -  *  fOii+Mi  +  O  +  B  (A19) 

o 

where  A  and  B  are  constants  analogous  to  integration  constants.  They  are  determined  by  the 
conditions 

%  =  0  (A20) 

and 

qm  =  0  (A21) 

from  which 


A  =  0 


(A22) 


and 


B  =  -  £i  IOn-i)0*i+^+l) 


(A23) 


Thus  after  simplification,  we  find 


% 


=  _  hi 

2m 


m 


m 


(m-n)  I  (2i-l)fc  +  n  I  (2m-2i  +  l)p. 
1  n  +  1 


(A24) 


From  the  nodal  displacements,  the  deflection  at  the  center  of  the  nth  element  is  calculated  to  be 


(A25) 


The  computer  algorithm  proceeds  as  follows: 

1.  The  strength  of  the  member  is  the  strength  of  the  weakest  element. 

2.  At  failing  load,  the  center  deflection  of  the  member  is  computed  from  (A23),  using  (A24)  and 
(A13). 


U.S.  Forest  Products  Laboratory 

Strength  of  Lumber  Under  Combined  Bending 
and  Compression,  by  John  J.  Zahn,  Res.  Pap.  FPL 
391,  FPL,  For.  Serv. ,  USDA.  13  p.  Madison,  Wis. 

Extensive  tests  of  western  2  by  6's  under 
eccentric  axial  load  support  the  conjecture  that 
adding  a  small  compressive  force  increases  bending 
strength.  Length  effect  was  tested  on  short  test 
members  by  simulating  long  members  with  a  finite 
element  computer  program.  A  small  sample  of 
long  members  was  tested  for  verification  of  the 
model  which  accurately  reproduces  the  mean  but 
underestimates  the  variance.  Further  work  on  the 
model  will  be  applicable  to  column  and  wood  truss 
design. 
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